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I derive expressions for various correlators of the gauge field and find propagators in a new gauge 
= 0. This gauge is a part of the wedge form of relativistic dynamics suggested earlier as the 
. . . tool for the study of quantum dynamics in collisions of hadrons at extremely high energies and in 

ultrarelativistic heavy ion collisions. The new gauge puts the quark and gluon fields of the colliding 
hadrons in one Hilbert space and thus allows one to avoid factorization. 

0\ . 

t^. I. INTRODUCTION 

^ : 

' In the previous paper |^ I explained physical motivation of the "wedge form of dynamics" as a promising tool to 
i explore the processes which take place during the collision of two heavy ions. In compliance with the general definition 
of dynamics given by Dirac |^ , the wedge form includes its specific definition of the quantum mechanical observables 
' on the space-like surfaces, as well as the means to describe evolution of the observables from the "earlier" space-like 
surface to the "later" one. Unlike the other forms, the wedge form explicitly refers to the two main geometrical 
I features of the phenomenon, strong localization of the initial interaction and, as a consequence, the absence of the 
. translational invariance in the temporal and longitudinal directions. 

Usually, every hadron (or nucleus) before the collision are considered separately, in their own infinite momentum 
frame. Thus, we always deal with two different Hamiltonian dynamics with their own definition of the time variable. 
This problem is seemingly resolved by the factorization scheme which replaces the true bound state by an artificial 
OS flux of free partons. However, the class of the physical processes which comply with this method is restricted to 
the inclusive production of jets with the high transverse momenta. Factorization is obviously violated already in the 
Qh] diffractive deep inelastic electron-proton scattering and, in fact, is not the case in any semi-inclusive process. 
JL . The constraints imposed by the factorization become critical when one is compelled to consider interaction of the 
Qj '' two bounded systems without appealing to the parton picture |Q, which is the case of the heavy-ion collision. One 
cannot design this experiment in such a way that the scale of the "hard probe," like, e.g., the dilepton mass in the 
Drell-Yan process, will become a measured parameter. In fact, the nuclei a priori probe each other at all scales and the 
expected factorization scale turns into ill-defined infra-red cut-off. Therefore, it is necessary to find a way to describe 
> both colliding systems, the hadrons or the nuclei, using the same Hamiltonian dynamic. This requirement follows 
?H ' solely from the fact that the definition of the field states (particles) depends on how the observables are defined. The 
problem is most painful for gluons since the choice of the gauge is one of the elements of the Hamiltonian dynamic. It 
manifestly affects the definition of the physical states of the gauge fields. Indeed, since every dynamics has a specific 
definition of the time variable. Gauss' law which defines the non-dynamical longitudinal fields, also looks differently. 

Factorization succeeds to fragment the cross section of the hadronic process into the hard partonic cross section, 
which is gauge invariant because the partons are considered as the free on mass shell particles, and the gauge invariant 
structure functions. By the failure of the factorization, one is lead back to the generic field-theory approach. It strongly 
requires the gluon fields from the incoming hadrons be given within the same gauge. 

The gluons are expected to be most abundant and active at the early stage of the heavy-ion collision. According to 
estimates of Shuryak and Xiong ||] most of the entropy should be produced due to the gluon interactions. However, 
before the collision, the quark and the gluon fields are assembled into the two coherent wave packets, the nuclei, and, 
therefore, the initial entropy equals zero. The coherence is lost and the entropy is created only due to interaction. 
Though it is tempting to rely on the general formula, S = Spplnp, which invariantly expresses the entropy via the 
density matrix p, at least one basis of states should be given explicitly. It is imperative to find what are the states 
that may be used for computing the entropy. 

In the wedge form of dynamics, the states of the quark and gluon fields are defined on the space-like hypersurfaces 
of the constant proper time r, = — z^. The states of scalar and fermion fields were discussed in Ref. In this 
paper, we continue the study for the gluons and augment our previous consideration by the gauge condition = 0. 
This simple idea solves several problems. First, it becomes possible to treat two different light-front gauges which 
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describe gluons from each nucleus of the initial state separately, as the two limits of this single gauge. Therefore, 
the new approach keeps important connections with the theory of the deep inelastic e-p scattering. This fact is 
vital for the subsequent calculations since the e-p deep inelastic scattering is the only existing source of the data on 
the nucleon structure in the high-energy collisions. [An alternative point of view is based on the classical model of 
the large nucleus in the infinite momentum frame M.] Second, after collision, this kind of a gauge becomes a local 
temporal axial gauge , thus providing a smooth transition to the Bjorken regime of the boost-invariant expansion. 

Most part of this paper is technica l, an d any relevant physical discussion of the results always goes only after their 
mathematical derivation. In section [I A, I derive equations of motion for the gauge field in the gauge = 0, find 
Hamiltonian variables and normaliza tion condition. Equ ations of motion are linearized and the modes of the free 
radiation field are obtained in section [IB . In section lIC the retarded propagator of the perturbation theory is found 
as the response function of the field on the external current. I carefully examine possibility to separate dynamics of 
the transverse and longitudinal fields and arrive at the negative conclusion. This part of calculation turned out to 
be the most time-consuming, since the gauge condition is inhomogeneous and none of the currently used methods 
is effective. The old-fasioned variation of parameters does work. To reass ure it s effectiveness, propagators of more 
familiar gauges, — and — 0, were computed in Ref. jo). In section II E I show that the previously obtained 
propagator solves the initial data problem for the gauge field. 

Section [II is devoted to the quantization of the vector field in the gauge A'^ — 0. I begin in [II A with computation 



of the Wightman functions and study the causal properties of the commutators in [II B . The latter appears to be 
abnormal, the Riemann function is not symmetric and penetrate the exterior of the light cone. However, behavior of 
the observables is fully causal and procedure of the canonical quantization is accomplished in [II C . Even though it 
is impossible to introduce transverse and longitudinal currents and thus fully separate dynamics of the corresponding 
fields, I found useful to classify various field patterns by the type of their propagation. Propagator of the transverse 
field is sensitive to the light cone boundaries while longitudinal and instantaneous parts of the field do not really 
propagate. These two fragments of the response function are derived in section ^ In section IV I study the limit 
behavior of the propagator in the central rapidity region and in the vicinity of the null-planes. I show that propagators 
of the gauges = and A* = 0, respectively are recovered. This result is practically important because it establishes 
connection of the new approach with the existing theory of the deeply inelastic processes at high energies. 



II. THE CLASSICAL TREATMENT 



A. Classical equations of motion 

Here we consider the case of pure glue-dynamics. We denote Afj,{x) — t"'A'^{x), the gluon field in the fundamental 
representation of the color group. Consequently, we have the field tensor, 

F^, = ^F;, = V^A, - V,A^ = df,A, ~ d,A^ - ig[A^, A,], 

where = 9^ — ig[A^{x)^ ...] is the covariant derivative on the local color group. The gauge invariant action of the 
theory looks as follows, 

5 = 1 Cix)d^x = J [-^g^\x)g''^ix)F,4x)FxA^) - rA^]y^d'x. (2.1) 

Its variation with respect to the gluon field yields the Lagrangian equations of motion, 

d,[{-g)'^h''Y^F,^] - tg{-g)^/'[A,,g'^Y^F,,] = {^gf/^f , (2.2) 

where is the color current of the fermion fields and g = detjg^^j. The equations are twice covariant: with 
respect to the gauge transformations in color space and the arbitrary transformations of the coordinates. In what 
follows we shall employ the special coordinates associated with the constant proper time hypersurfaces inside the light 
cone of the collision point t = z — Q. The new coordinates parameterize the Minkowsky coordinates {t,x,y,z) as 
(rcosh?7,a;,y,rsinh?7). In addition, we impose the gauge condition Ar — 0. The corresponding gauge transformation 
is well defined. Indeed, let A^{x) be an arbitrary field configuration and A'^{x) its gauge transform with the generator 

UiT,r],rl)=PrC^v{- f Ar{T',r^,rl)dT'} (2.3) 
Jo 

Then the new field, A'^ — UAfJJ^^ + d^UU^^, obey the condition A'^ — 0. Imposing this gauge condition we arrive 
at the system of four equations: 
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C{x) ^ -dndrA,^ + rdrdrAr - ig{-[A,^,drA,j] + T[Ar,drAr]} - rf ^ 0, (2.4) 
T r 

^ drTdrAr + —djj{d,^Ar ~ BrA,^) + TdsipsA^ — drAg) 
T 

-ig{-d„[A^,Ar] + Tds[As,Ar] + -[Arj, F,,r] + t[A,, F,,]} - Tf = , (2.5) 
T T 

- dr-drA„ + -dr{drA„ - d„Ar) ~ ig\-dr\Ar,AA + -\Ar, FrS\ ~ rf = (2.6) 

r T T T 

Here, we use the latin indices from r to w for transverse x- and y-components, r, w = 1,2. We shall also use the 
arrows over the letters to denote the two-dimensional vectors, like k = {kx,ky), \k\ — k±. The latin indices from i 
to n, i,...,n = 1,2,3, will be used for the three-dimensional internal coordinates = {x,y,i]) on the hypersurface 
T = const. The metric tensor has only diagonal components grr — —§,xx — ^Syy — 1? Svn ~ ^^'^ &rst of 

these equations contains no second order time derivatives and is a constraint rather than a dynamical equation. The 
constraint weakly equals to zero in classical Hamiltonian dynamics and serves as the condition on physical states in 
quantum theory. The canonical momenta of the theory are as follows: 

1 1 • • 

=0, = -Fr^ = - A„ Pr = TFrr ^ T Ar ■ (2.7) 

T T 

Hereon, the dot above the letter denotes derivative with respect to the Hamiltonian time t. Because of the gauge con- 
dition, the canonical momenta do not contain the color commutators. After excluding the velocities, the Hamiltonian 
can be written down in the canonical variables: 

H = J dl^drlT{]^pyi + + ^F^rF^r + \FrsFrs + f A^ + f A,.} (2.8) 



Then the equations (2.5) and (2.6) are immediately recognized as the Hamiltonian equations of motion. The Poison 
bracket of the constraint C with the Hamiltonian vanishes thus creating the generator of the residual gauge transfor- 
mations which are tangent to the hypersurface. Conservation of the constraint is a direct consequence of the Lagrange 
(or Hamiltonian) classical equations of motion as well. 

The normalization condition for the one-particle solutions is obviously derived from the charge conservation law. 
For the gauge field this is impossible. Therefore we shall accept the condition which supports self-adjointness of the 
homogeneous system after its linearization. This leads to a natural definition for the scalar product of the states of 
the vector field in the gauge = : 

{V,W)^ f dij [ d?rTg'''V*idrWk (2.9) 



where g^^ is the metric tensor of the three-dimensional internal geometry of the hypersurface r = const . This norm 
of the one-particle states prevents them from the flow out of the interior of the past and future light wedges of the 
interaction plane. 



B. Modes of the free radiation field. Field of the static source 



As a tool for the future development of the perturbation theory, we need to find the propagators and Wightman 
functions when the nonlinear self-interaction of the gluon field is switched off. In this case the system of the equations 
for the nonvanishing components of the vector potential and the constraint look as follows: 

[drTdr - -dl - ra2]A, + d^rd^A, + -d,,A^\ = -Tf , (2.10) 

T T 

\d^\d^-\df\A^^-d^d,A,^-T3'' , (2.11) 
r r T 
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C{x) = -drjdrArj + rdrOrAr - rf = . (2.12) 
T 

where includes all kinds of the color currents. An explicit form of the solution for the homogeneous system is found 
in the Appendix 1. In compliance with the gauge condition which explicitly eliminates one of four field components 
we find three modes V'^'^^ of the free vector field. Two transverse modes obey Gauss law without the charge and have 



the unit norm (see Appendix 1) with respect to the scalar product (2.9) 



p-TTv/2 I ky \ ^ p--^vl2 



'^kyR^X-^.ik±r) (2.13) 



The mode y'-^-' is constructed from the functions R^J''_j_^{k±T) = R^^^_^i^{k±T\s) corresponding to the boundary 
condition of vanishing gauge field at t = 0. This guarantee continuous behavior of the field at r = 0. Indeed, at 
T ^ 0, the normal and the tangent directions become degenerate. As long as A'^ = is the gauge condition, the 
continuity requires that A'' ^ at r ^ 0. 

In order to simplify some subsequent calculations it is useful to write down explicitly the physical components of 

electric and magnetic fields of these modes, 5™ = a/— gg™" A71 and = — 2^^e™'"F/„; 

H^_^l{k^T)e"''>+'''^ . (2.14) 

Therefore, the mode V*^^^ can be obtained from the mode V'^^'' by simple interchange of its electric and magnetic 
fields. Using the standard wave-guide terminology, one may call mode V*^^^ as the "transverse electric mode" and the 
mode T^^^-* - as the "transverse magnetic mode" . 

An equivalent full set of the transverse modes carries instead of the boost v, the quantum number 6, rapidity: 
^0 = k± cosh0, fcs = k±_ sinh0. These functions can be obtained by means of the Fourier transform, 

and have the following form, 

- ( ) e--— .g(x) = j-J— I fc,/r | e^^ , (2.16) 





where 



fi{T,v) = k± smh{0 -V) f e-''=^^'™''^(''~'')dr' = ztanh(6' - ^)(e-^*^-L^™>'ii(e-'7) _ l) 
Jo 



^—ik±T cosh{9 — ri) ^~ik±T cosh{8~rj) 



Mr.,) = kl I e-.— = ^'"^r ^^^^^^^^ . (2.17) 



The norm of the Coulomb mode V^^\ as defined by Eq. ( |2.9| ), equals zero, and it is orthogonal to V^^'' and V^'^K 
Though this solution obeys equations of motion without the current, it does not obey the Gauss law without a charge. 
Therefore, it should be discarded in decomposition of the radiation field. However, it should have been kept if we 
considered the radiation field in the presence of the static source with the r-independent density p{k, u) = rjZ (r) = 
const{T). In this case its definition can be completed using the Gauss law: 

y(3)(^) ^ ( krQ-iMk^r) \ 
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The coordinate form of this solution is noteworthy. The physical components, E™' = y/—gg™''Ai , of the electric field 
of the "r-static" source can be written down in the integral form, 



with the kernel 



^lstat)i'^^^i^Vi) = / dr2d-q2Ki{T\ri -^2,771 - m) pir2 , V2) , (2.19) 



where 



R12 = (r-i + sinh^ ,y)l/2 = [(^- _ )2 _ ^2 ^^^^^^2^^^ _ ^^)ji/2 ^ ^2.21) 



is the distance between the points ri,7yi and 7^2,772 in the internal geometry of the surface r — const. The technical 
details of derivation of the last expression will be adduced in Sec. ^ Eq. ( |2.20| ) is an analogue of the Coulomb law 



of electrostatics, except that now the source has a density which is static with respect to the Hamiltonian time r. In 
fact, the source is static if it expands in such a way that Tj^(r, 77, r) does not depend on r. These expressions will be 
helpful in recognizing the origin of various terms in the full propagator which is calculated below. 



C. Propagator in the gauge A'^ = 



Calculation of the propagator in the gauge A'^ = associated with the system of the curved surfaces r = const 
meets several problems. Three methods are commonly used in the field theory. One of them strongly appeals to the 
Fourier analysis in the plane Minkowsky space which is not applicable now because the metric itself is coordinate- 
dependent. The second method uses the path- integral formulation which is also ineffective because of the explicit 
coordinate dependence of the gauge-fixing term in the Lagrangian. One could also try to study the spectrum of the 
matrix differential operator, to find its eigcn-functions and to use the standard expression for the resolvent. However, 
the extension of the system for the non-zero eigenvalues leads to the unwieldy equations. On the other hand, the Green 
function of the perturbation theory must coincide with the one which solves the problem of the gauge field interaction 
with the classical "external" current. For this reason, we shall compute the Green function in a most straightforward 
way: we shall look for the partial solution of the inhomogeneous system using the method of "variation of parameters" . 

Let us start derivation of the propagator in the gauge A'^ = with obtaining the separate differential equations 

for the 7y-component of the magnetic field, 5* = dyA^ — dxAy-, the transverse divergence of the electric field, (p = 

• • • 

T{dxAx + 9yAy)', and the ry-component of the electric field, a ~ Ati/t . In terms of the Fourier components with 

respect to the spacial coordinates, these equations read as 

{dr + \dr+'^ + fci]*,v(r) = 7., r) , (2.22) 



[drTdr + —]f{k,iy,T) - iTiykla{k,i^,T) = -dr[T^j^{k,iy,T)] , (2.23) 
r 

[drTdr +Tkl]a{k,iy,T)~—^(k,u,T) = -dr[T^f{k,iy,T)] , (2.24) 

T 

where j'^ = dyjx — dxjy, j"^ — dxjx + dyjy. Using the constraint conservation, which may be explicitly integrated to 

^p{k,v,T)) + iva{k,v,T) — T]'^ [k,v,T) ^ — PQ{k,v) = const(T) , (2.25) 
one easily obtains two independent equations for (p{k, r) and a(fc, r): 

[dl + -dr+'^ + kl\^{k, V, t) = kl[p{k, r) - po(fc, v)\ ~ -drir^j^k, r)) ee ^ , (2.26) 
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[a^ + la, + 4 + kl]a{k, V, t) - -^[p{k, V, t) - po{k, ly)] - -dr{T^f{k, v, r)) = /" , (2.27) 

T T- T T 

The constant of integration po(^i has meaning of the arbitrary static charge density and it should be retained until 
the Gauss law is explicitly imposed on the solution. In what follows, we shall not write it down explicitly assuming 
that it is included in the true charge density p(fc, j/, r). Equations ( ^.22 ), ( ^.26| ) and ( 2.27 ) can be solved by the 
method of "variation of parameters" : 

^(r) = - T2dT2n{T,T2)f{T2) , (2.28) 

where T stands for anyone of the unknown functions in these equations, and / for the corresponding right hand side. 
The kernel 

H(r, r2) = H'i) ik^T)Hl^^ (fc^r^) - H^^ ik^T)H^/J [k^r^] . 

is a usual bilinear form built from the linearly independent solutions of the homogeneous equation. [ The Wronskian 
of these solutions is exactly 4/i7rT2. ] Taking = we obtain the first equation for the components A^lk, v, t) and 
Ay{k, t) of the vector potential: 

*(fc,z/,ri) = i[-kyA, + k.,Ay] = — / T2dr2H(Ti,T2) I [-kyf{T2) + k,jy{T2)] . (2.29) 

4 Jo 

In order to find the second equation for the x- and y-components and the equation for A^(fc, r) we must integrate 
twice: 

in f^^ dr' r' ^ 1 ^ 

^k,y,n) = i[k^A^+kyAy]^ — — 'H(T',T2)T2dT2[-kip{k,U,T2) + -dr,{Tlf{k,V,T2))], (2.30) 

4 Jo '''Jo 



A^{k, V, Ti) = — / r'dr' / H(t', r2)r2dT2 [— p(fc, T2) + -9., (r|j''(fc, T2))\ . (2.31) 

4 Jo Jo '^2 '''2 

Here, the integration over T2 recovers the electric fields at the moment r', whilst the integration over t' gives the 
vector potential at the moment ti . It is convenient to start with the latter one which has the limits T2 < t' < ti. Let 
us consider t he ma in line of calculations in detail, using the 77-component as an example. The first integration follows 
the formula ( |A2.l|) , 

fc'^+i r\rrH\^{k^T')dr' = Ri^lik^r^ - R%{k^r2) , (2.32) 



and the terms emerging from the lower limit T2 can be conveniently transformed according to the relation (see 
Appendix 2), 

R%{k^T2)HZ\k^T2) - R%{k^T2)H'll\k^T2) - ^S^...(fc±T2) (2.33) 

As a result, one obtains, e.^., the following formula for A^(fc, j^, ri): 

= ^ r T2dT2[R[%ik^T,)Hl^\k^T2) - ^'l (fc^n )i? f^^^ (fc^r2 ) - 1 Si,,, (A;^r2 )] 

x[^p{k,,.,T2) + -dr,{T^f{k,,.,r2))] . (2.34) 
T2 T2 

In order to eliminate the charge density p from the integrand and to separate the transverse and the longitudinal parts 
of the propagator, all the terms of this formula should be integrated by parts with explicit account for the charge 
conservation which reads as 

lT[k^f{k, V, t) + k^f{k, V, t) + v3-^{k, V, t)\ + drp{k, z., t) = . (2.35) 
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We have in sequence: 



IV I p(T2)n-^^' [k±T2) ^ IV ; P[T2)aT2 

lo "^2 Jo aT2 

Jo 



(2.36) 



zvj^ dT,H^(k^T^d,,{^]^{T^)=Tlf{T,)H^(k^T,)+ / T^dr^lR^flik ^T^) + v^ R%^{k ^T^)]^ (t^) . (2.37) 
In a similar way we have, 



p{T2)dT2 



IV I p{T2)si.tu{kl_T2) ^ iv . 

lo ^2 Jo dT2 

ivQ-l,iu{k^Tl)p{Tl) -V I T2dT2Q-l.ty{k^T2)[kxf{T2) + kyf{T2) + vf{T2)] 

Jo 



(2.38) 



' dT2S,,Uk±T2)dr,{T^f{T2)) - T^'" (n )si,, (fc^n ) + V^ f ^ TadTa [g_l,,,(fciT2) - Ql,.. (fc^Tz )] j" (t2 ) . (2.39) 

Jo 

Assembling these pieces together and repeating the same calculations for the function <i> one obtains three different 
terms which contribute to the field A produced by the current j: A = A^*^'' + A*^^^ + ^^™^*). 



The transverse field A^*''^ is defined by the integral terms in the R.H.S. of Eqs. (2.36) and (2.37). It can be 
conveniently written down in the following form: 



A[''-\x,)^ / d'x2e{T^-T2)A'^\xi,X2)r{x2) 



where 



^ — OO \ 1 o 



(2.40) 



(2.41) 



A=l,2 



can be easily recognized as the Riemann function of the original homogeneous hyperbolic system. The Riemann 
function solves the boundary value problem for the evolution of the f ree r adiation field. It is obtained immediately 
in the form of the bilinear expansion over the full set of the solutions ( 2.13 ) of the homogeneous system. In fact, this 
is a sole evidence that A^*'") may be associated with the transverse part of the propagator. Then the rest part is the 

propagator (response function) for the longitudinal field. 

The dynamical longitudinal field A^-^^ originates from the integral terms in the R.H.S. of Eqs. ( 2.3§| ) and (2.3£): 



M^Hn^m^n) = J T2dT2 J dri2d^f2A^i!^{T2;'ni - T]2,ri - r2)j"(T2, ??2, 7^2) 



(2.42) 



and the kernel of this representation. 



A[^^^(r2;r7i -772,^1-^2) 



r dvd^k 






ksQ~l,iu{k±T2) 


^iv{r]i-r)2)+ik[ri-r2) 


1 {27Trkl 


V 


I 


vQlSu{k±T2) 


rn 



(2.43) 



does not allow for the bilinear expan sion w ith two temporal arguments, and, as we shall see in a while, the retarded 
character of the integration in Eq. (2.42) is not sensitive to the light cone boundaries. In fact, the electric field 

(L) . 

El — Ai is simultaneous with the current j 



The last one, instantaneous, part of the solution comes from the boundary terms in Eqs. (2.36)^( 2.391 ) which were 
generated via integration by parts. It depends on a single time variable ri. Using two functional relations, (2.33) and 



ITT V- 



dx 



- — [Qi,iiy{x) - Q~i,i„{x)] 

ITT 



(2.44) 
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( see Appendix 2 ) its Fourier transform can be presented in the form 



krQ-lAu{k±Ti) 
yQl.:iv{k±Ti) 



(27r)3ifc2 

which leads to the Poisson-type integral, 

^m"''*'' (n, ^i) = / dr2diq2K.m{Ti\ri ~ ^2, 771 - ?72)/0(ti, ^2, 772) , 



with the instantaneous kernel, 

A^m(T;r, 77) 



dvdk 6"^''+*'=^"* 
(27r)3 ik\ 



krQ-l.iv{kA.T) 
vQl^tu{k±T) 



(2.45) 



(2.46) 



(2.47) 



Therefore, this term represents instantaneous distribution of the potential at the moment ri, corresponding to the 



charge density taken at the same moment. Recalling that the charge density p{k,v,Ti) in Eq. ( 2.45| ) still includes 
the arbitrary constant po{k^ v), we see that imposing the constraint indeed affects only the potential of static charge 
distribution and puts it in agreement with the Gauss law. If po is set to zero, then the gauge is completely fixed and 
all calculations with this propagator will be gauge invariant. 

For the practical calculation, it is easier to work with the components of the current rather than to keep the charge 
density p in its original form. Otherwise, e.g., the expression for the vertex function will be unwieldy. Replacement 
follows the prescription: 



p(Ti,?y2,r2) 



dp 



r2)] 



and restores an extraneous "initial" configuration of the static charge which has been previously removed from 
Eq. ( ^.45| ). Being r-dependent, its vector potential cannot be eliminated by the residual gauge transformation. 
However, this artificial contribution corresponds to the easily recognizable static pattern in the longitudinal part of 
the propagator and is under the full control. Keeping this fragment in mind, we arrive at the standard form of the 
representation, 



A'™'**^(Tl,?7i,fi) = j T2dT2 y"d?72dV2A|™'**^(Ti;?7i - 772, fi - f2)j'"(T2,772,f2) 

with the kernel given by the formula. 



(inst) 



{ti;Vi - V2,ri - ^2) 



diycPk 
(27r)3fc2 



krQ-lauik±Ti) 
vQl.iu{kl_Ti) 





■fc/ 


I 


V 



^iv{rii-ri2)+ik(ri-r2) 



(2.49) 



Eqs. ( 2.40| )-( 2.43 ) and ( 2.48 ), ( 2.49 ) present the propagator in a split form. Different constituents of this form 
are preliminary identified as transverse, longitudinal and instantaneous parts of the propagator. It would be useful 
to learn if the same kind of splitting is possible for the current itself. An affirmative answer (as in the cases of the 
Coulomb and radiation gauges) would be helpful for the design of the perturbation theory. To answer the question, 
one should substitute the different pieces of the solution into the left hand side of the original system of differential 
equations. This leads to the following expressions for the Fourier components of the three currents: 



Tj7*V)(fc,^^;r)=r/ 



(fc, j^;t) + -V 



kr^l.iu {k±T^ 

vk\s-i^iu{k±_T) 



dp 



V d 



krT^P 



V 



T{k^f + kyjy) 

[Q-i.,Ak±T){k^f + kyf) + Qi,Uk±r)f] 



(2.50) 



(2.51) 



-1 d 
ik"] dr 



krTQ-iai,{k^T) 
VT^^Ql,iu{k^T) 



rdp\ 


1 


dr J 





krSl,i„ik±T) 

vk\s^i^iy{k^T) 



dp 



(2.52) 
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Providing the current is conserved, these three currents, added together, give the full current from the right hand side 
of the system. Therefore, the solution is correct. However, none of these three currents carries any signature of being 
longitudinal or transversal in a usual sense. None of them has zero divergence, since the operator of the divergence 
does not commute with the differential operator of the system. No desired simplification is possible in our case. 

In fact, even the above splitting of the potential has no real physical meaning. To see it explicitly, let us find the 
divergence of the electric field, divE — dm£™ [again, for brevity, in the Fourier representation]: 

divE'-'^Hk, v-t) = i[Q^iMk±T) - QiMk±r)]{:yT^f - ^{k^f + kyjV)) , (2.53) 
divE(^)(fc,z.;T) =z(r2 + ^)[(A:,r + fc,j^')Q_i,,,(fc^T) + z.j''Qi,,,(fcir)] , (2.54) 
divE(™^*)(^,i.;T) =p(fc,z.;T)-z(T2Q_i,,,(A:ir)-^Qi,,,(fc^r))[(A;,r + fcyj") + i^j''] • (2.55) 



Only the divergence of the true retarded component of the field E*^*''-' turns out to be zero. The term which prevents 
divE^*''^ from being zero is due to non-symmetry of the propagator. A'''' ^ A'"''. It appears when the 0- function in 



Eq. (2.40) is differentiated with respect to Hamiltonian time r. This term is vital for obtaining the expression that 
obeys Gauss law constraint, divE(fc, v] r) = p{k^ i^',t). [ We remind that p(fc, ly; r) still includes an arbitrary constant 
Po{k,i^).] 

The known examples, when the transverse and the longitudinal fields are separated at the level of equations of 
motion, are related to the narrow class of homogeneous gauges. Impossibility to perform a universal separation of the 
transverse and longitudinal fields thus appears to be a rule rather than exception. It reflects a general principle: the 
radiation field created at some time interval has the preceeding and the subsequent configurations of the longitudinal 
field as the boundary condition. Dynamics of the longitudinal field falls out of any scattering problem in its S'-matrix 
formulation. However, this dynamics is, in fact, a subject of the QCD evolution in deep inelastic scattering 

D. Initial data problem in the gauge = 

We obtained the expression for the (retarded) propagator as the response function between the "external" current 
and the potential of the gauge field. We must verify that the same propagator solves the Cauchy problem for the 
gauge field. This can be easily done by presenting the initial data at the surface t = tq in the form of the source 
density at the hypersurface r = tq. 



V^J^T2) = v/^i(^g""(ro)[5'(r2 - ro)A™(f, r?) + S{t2 - To)A:„(f, ^)\ , (2.56) 
where A,„(r, 77) and A'^^r, rj) are the initial data for the potential and its normal derivative on the hypersurface r — tq 



Usua Uy, i t is assumed that the real currents vanish at r < tq. Substituting this source into the Eqs. (2.4C), (2.42) 



and ( 2.46 ), and taking the limit of ti — > tq, we may verify that the standard prescription for the solution of the initial 



data problem, 

Mxi)= f dV2 d772 Ah„(a;i,a;2)T^A'"(.T2) , (2.57) 

holds with the same propagator lS.im{xi^X2) that was used to solve the emission problem. For example, in the limit 
of r ^ To, the ?7-component of the vector potential is a sum of three terms. 



?l{k^ro)H[l\k^To) - R^^l{k^r,)H\l\ 
roiy[R'flik^To)R%,{k^ro) - R[%{k^ro)R%,Jk^^^^^ } , (2.58) 



A^v'-Hro + 0) = ^ { [R[%{k±ro)Hl'J{k^ro) - R[]Uk^ro)Hl'J {k^To)][,.A^ - klA,] 



AII^Hto + O)^ -si.,^(fciTo)A0 + roQ_i.„,(fc_LTo)A;-^s_i,,^(fc_Lro)^^ + — Qi,,^(/ciro)^'^ } , (2.59) 

Kj_ Tq Tq 
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4"^*)(ro + 0) = ^Qi^„,(fc^ro)[To^; 



To 



(2.60) 



where we have denoted: Aa 



^y^^y Here, Eq. ( 2.6C ) follows from Eq. ( 2.46 ) and takes care o f the consis tency 
between the charge d ensit y at the moment tq and the initial data for the gauge field. Using relations (2.33) and (2.45|) 



and adding up Eqs. (|2!5^)- (^^ 

we come to a desired identity, A,,(ro + 0) = 



When the initial data Am{f, rf) and A'^{r, rj) cor respo nd to the free radiation field, then only the part A|^^(a:i, 2:2) 
of the full propagator "works" here, and only Eq. (2.58) may be retained. The other two equations acquire status of 
the c onstr aints imposed on the initial data. Since the current is absent, we have A^^^^ = on the left hand side of the 
Eqs. (2.5£). Then the right hand side confirms that the kernel JC is orthogonal to the free radiation field modes. Since 
the charge density p vanishes, we have ^("'**) = 0, which is equivalent to the Gauss law for the free gauge field. The 
two transverse modes already obey these constraints. This fact provides a reliable footing for canonical quantization 
of the free field in the gauge A"^ — 0. Indeed, the Riemann function induces commutation relations for the gauge 
field. It can be found via its bilinear decomposition over the physical modes. Thus one can avoid technical problems 
of inverting the constraint equations. (See Sec. III.) The longitudinal part of the propagator will be studied in details 
in Sec. fvl. 



E. Gluon vertices in the gauge A'^ — 

The terms proportional to the first and the second powers of the coupling constant in the classical wave equations 
may be viewed as the external current and allow one to define the explicit form of the 3- and 4-gluon vertices. One 
should start from the solution of the Maxwell equations, 

Ai{z,) = j d^x^i,l{z^,x)./^)J^{x) . (2.61) 

with the color current of the form 

v/^gMJa' (2;) - -.g/afccy^^g'="(x)g™'(x)[a,„(Af(x)A^(a;)) + Af(x)5„A,=,(a;) + A:;(x)9„A^(a;)] 

-.9'v/^i(^/afcc/c<i/.g'="(x)g'"'(x)Af(a;)Al(a;)A:j(x) . (2.62) 

In perturbation calculations, every field A{x) in the RHS of this expression is a part of some correlator A(x, zjv). 
The components of the metric depend only on the time r while the derivatives affect only the spacial directions 
m" = (r, rj). Moreover, in these directions, all the gluon correlators depend only on the differences of the coordinates 
and can be rewritten in terms of their spacial Fourier components. After symmetrization over the outer arguments 
zisfi one immediately obtains, 

Kt"(Pi,P2,P3; r) = -irfabc 5{p, +P2+ P3) [g'"(P2 - Ps)" + g"'(p3 - Pi)' + g"' {pi ^ P2)"] , (2.63) 

where p" — g"'^Pk, and the components of the momentum in the curvilinear coordinates are equal to pk — {px,Py, v). 
The four-gluon vertex has no derivatives and is the same as usually. 



III. QUANTIZATION 

The second quantization of the field has several practical goals. We would like to have an expansion of the operator 
of the free gluon field like 

A,{x)^ J2 /d^fcd^[cA(z^,fc)T/Jg(a;)+4(j.,fc)V;(^^);(x)] , (3.1) 

with the creation and annihilation operators which obey the commutation relations 

[cx{'y,k),cl,{iy',k')]=Sxy6{iy-iy')S{k-k'), [cx{iy,k),cy{iy',k')] = [c{{iy,k),c{,{iy',k')]^Q . (3.2) 

Once obtained, commutation relations ( |3.2| ) allow one to find various correlators of the free gluon field as the averages 
of the binary operator products over the state of the perturbative vacuum and express them via the solutions V^^\{x) 

('2) 

and V J (x). For example, the Wiehtman functions, 

h>k\i 
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iAlo,^J{x,y) = {0\A,{x)Aj{y)\0) = 



A=l,2 



(3.3) 



serve as the projectors onto the space of the on-mass-sheh gluons and should be known exphcitly in order to have a 
good definition for the production rate of the gluons in the final states. With these two Wightman functions at hand, 
one immediately obtains the expression for the commutator of the free field operators, 



^Q,tj{x,y) = -i{Q\[A,{x),Aj{y)]\Q) = Aio,,j{x,y) ~ Aoi,,j(a;,y) 



(3.4) 



which should coincide with the Riemann function of the homogeneous field equations. The program of secondary 
quantization do not rev eal a ny technical problems if we give preference to the holomorphic quantization which starts 
commutation relations (B^) for the Fock operators. However, if we prefer to start with the canonical commutation 
relations for the field coordinates and momenta, then one should postulate them and derive ( |3.2| ) as the consequence. 

The way to obtain the canonical commutation relations in cases of the scalar and the spinor fields is quite straight- 
forward. For the vector gauge field we meet a well known problem, an excess of the number of the components of the 
vector field over the number of the physical degrees of freedom. For example, in the so-called radiation gauge, A^ = 
and divA = 0, we write the canonical commutation relations in the following form M, 



[A{^,t),E,{y,t)]^6l^{^-y)= J 



(2^)3 V ' k2 
[^,(x,t),^,(y,i)] = [E,i^,t),Ej{y,t)] =0 



(3.5) 



thus eliminating the longitudinally polarized photons from the dynamical degrees of freedom. The function plays a 
role as the unit operator in the space of the physical states. Here, i,j = 1, 2, 3 and the number of relations postulated 
by equations ( |3.5| ) apparently exceeds the actual number required by the count of the independent degrees of freedom, 
A = 1, 2, of the free gauge field. The Fourier transform of the function (5*J is easily guessed because the basis of the 
plane- wave solutions is very simple , and it can be obtained rigorously by solving the system of constraint equations 
0,0. A similar guess or procedure in our case is not so obvious. We have the gauge condition = as the 
primary constraint and the Gauss law as the secondary one. The latter can be resolved in a way which allows one 
to exclude the ry-components of the potential and the electric field from the set of independent canonical variables. 
Thus, only x- and y-components are subject for the canonical commutation relations. To resolve the constraints, one 
anyway needs the integral operators with the kernels built from the so lutions of the Maxwell equations in the gauge 
A'^ = 0. Therefore we shall proceed in two steps. In section IH A we shall sketch the results for the Wightman 
functions (3.2). These, will be used for the explicit calculation of the free field commutator (3.4) in section [HB and 
for the study of its causal behavior. 



A. Gluon correlators in the gauge A'^ — 

Here, we shall write down components of the field correlator Aio,ij {x, y) in the curvilinear coordinates u = (t, 77, r). 
We shall denote their covariant components as Aio.ifc(Mi, 1*2). Later we shall transform them to the standard 
Minkowsky coordinates and find the correlators of the temporal axial and the null-plane gauges as their limits in 
the central rapidity region and in the vicinity of the null-planes, respectively. The most convenient for this purpose 
basis consists of the transverse modes v'^^\ The mode v'^'^^ gives the following contribution to the correlator Aio^ifc: 

iA*-^^ (1 2) = / — I '^^^ ^'^"'^^^"^"^" ^^fc(rl-f2)^-^fcxTlCosh(6l~^;l)+^fc^r2Cosh(6)-^;2) ^3 g-j 

J-00 2 J (27r)3 

Realising that d9/2 = dk^/2k", we recognize a standard representation of this part of the correlator in terms of the 
on-mass-shell plane waves decomposition. 

The second part of the correlator is determined by the mode v'^^^ and has the following components: 

AfoU(l>2) ^^^^ ^1 ^^e^^>--)A(^,ri,r;i)/r(0,r2,ry2) , (3.7) 



— OO 



J (|^e^''''^"''^^^/i(^'^i'^i)/2*(^.^2,^2) = zAPol„.(2,l) , (3.8) 
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^ yy (2^)3 f2{e,Tum)m,T,,,n) ■ (3.9) 

One may easily see that all components of Aio(l, 2) vanish when either ri or T2 go to zero. 



B. Causal properties of the field commutators in the gauge = 



Ca usal properties of the radiation field commutator may be studied starting from the representation (p. 41) . Using 
Eqs. (3.6) and (3.7) we may conveniently wright contribution of two transverse modes in the following form: 



kl 



dO sin A; I $ 



(3.10) 



d k kip k 



(27r)3 kl 



de 



cosh 2rj 



sinh^ 6 — cosh^ rj 



(sinA:^$ - sinA:^$i + sinA:x$2) , (3-11) 



where we ha ve int roduc ed th e following notation: 2r] — rji — r]2 , r = r\ — , — t\ cosh(0 — jyi) , <i> = $1 — $2 . 
The sum of ( ^l0| ) and (^Hl) can be rearranged to the following form, 



iAo,rs(l, 2) = i 



d^k de 

Tl /•-r2 



ikr 



5r^ sin fc I $ 



kf kg 



[sinfc^<I>i — sinA:^$2] 



+krks cosh{rii — ri2) / dr' / dr" sin[fc_Lr' cosh(0 — ry) — A:^r' cosh(0 + 77)] 
"'0 Jo 



(3.12) 



Joining the integration d^k dO into the three dimensional integration (i"^k/|k| in the Cartesian coordinates, the first 



integral in ( 3.12 ) 



^0(1,2) 



d'^k de 



e'''''sin A;_L$ 



(2^) 



sign(ti - ^2) 
27r 



5[{h-hf-{v,-V2f] , 



(3.13) 



is easy to calculate, and to recognize it as the commutator of the massless scalar field. It differs from zero only if the 
line between the points xi and X2 has the light-like direction. In this way we integrate the first and the third terms 
in the integrand of the Eq. (3.12). To reduce two integrals in the second term to the same type, we must exclude the 
factor 1/fc^ using the fundamental solution of the two-dimensional Laplace operator. 



Kj_ 



d^ 

yr 



In 1^ — r^e 



(3.14) 



After that we arrive at the final result, 



Ao,.,(l, 2) = -drsDoil, 2) - cosh(r;i - r,2)drds f dn f rfT2i?o(l, 2) 

Jo Jo 



+drds 



(27r)2 



ln|e-r|[<5(r2-e')-5(r|-^ 



-2 



(3.15) 



From this form it immediately follows that the commutator of the potentials vanishes at ti = T2- Even more strong 
result takes place for the commutator of the two electric fields. 



[Er{l),E,{2)] = ,:^^zAo,..(l,2) = 

OTIOT2 



— 5r 



92 



dTldT2 



cosh(77i - 772) 



92 



zZ?o(l,2) 



(3.16) 



This commutator vanishes everywhere except for the light cone, in full compliance with the microcausality principle 
for the electric field which is an observable. However, this does not happens for the commutator of the potentials 
since they are defined nonlocally. It does not vanish neither at space-like nor at space-like separation, because the 
line of integration which recovers the potential at the point X2, in general, intersects {e.g. at some point X3) with the 
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light cone which has its vertex at the point xi, and the commutator of the electric fields at the points xi and ^3 is 
not zero. 

Similar results take place for the commutator of the r/-components of the potential and the electric field. The field 
commutator, 

[i?^(l),£;,(2)] = ^-^iAo,^,(l,2) = -iVii^o(l,2) , (3.17) 

is entirely causal, while the commutator of the potentials, 

[A^{1),A^{2)] = iAo,^^(l, 2) = -iVi / ' ridn / ' rsdrsi^oll, 2) , (3.18) 

Jo Jo 

does not vanishes at space- like distances, except for ti = T2. Finally, the formally designed commutator between the 
r- and 77-components of the electric field, the two observables, 

[EAl),E,i2)] = ^^Ao,.,(l,2) ^ -^(^^0(1,2)) , (3.19) 

is entirely confined to the light cone, while the commutator of the potentials, which are not the observables, 

K(1),A,(2)] =zAo^,,,(l,2) = -y^ — i r2dT2^^Do{l,2) , (3.20) 

does not vanish at the space-like distance, even at ti — T2. This result, however, is not a subject for any concern 
since the potentials are defined nonlocally and commutation relations for electric and magnetic [cf. ( |2.14| )] fields are 
reproduced correctly. Moreover, we have argued above that the T^-components of A and E are not the canonical 
variables since the constraints express them via x- and y-components. 

The "acausal" behavior of the Riemann function Aq''(1,2) may cause doubts if the gauge A'^ = allows for 
meaningful retarded and advanced Green functions which, by causality, should vanish at space-like distances. Fortu- 
nately, this anomalous behavior appears only for the gauge-variant potential; th e resp o nse fu nction s for observable 



electric and magnetic fields are causal. This can be easily seen, e.g., from Eqs. ( 2.22| ), ( 2.2t ) and ( 2.27 ), the usual 



inhomogeneous relativistic wave equations for various physical components of the field strengths £ and B. 

C. Canonical commutation relations in the gauge A'^ — 



A proof of the commutation relations (3.2) for the Fock operators follows the standard guidelines |9|. The creation 



and annihilation operators are defined via relations, 

/• •(■^)* 
d^^[V^^*{x)Ad^,T) - V.j:.,,{x)A{^,T)] , 

c1(^., k) ^ {A, V^^^) = zg'^- J rf=^x[A,(x, t)vX{x) - l(x, r)V^~l{x)] . (3.21) 

The latter result in the following expression for the commutator, 

dWyg''^(x)g'"(y) {[A,(x,T), A(y,T)] (y.j:,,{x)V^^Jy) - V^^p.^^xW^^liy)) 

,(A)* ,(A') , , , 

+ [A{^,T),Aiy,T)] n£;,(^)V^.'^';Jy) + U.(x,r),Az(y,r)] l^iV*(^)l^j,-'„(y)} ■ (3-22) 

Most of the terms in the second line vanish due to the commutation relations. Next, we rely on the following guess 
about the form of the commutator, 

[A,{x),AM]= E f (^ISS'^^^l^^^^^^^lSl^'^^^lS^^y^ ' (3.23) 

which leads to the proper equal-time commutation relations for the independent canonical variables. Finally, explicitly 



using the orthogonality relations for the eigenmodes V'^^\ we immediately obtain the commutation relations (3.2) 
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IV. LONGITUDINAL PROPAGATOR AND STATIC FIELDS 



In this section we shall find the explicit expressions for the kernels (2.43) and (2.47) which represent the longitudinal 
and instantaneous components of the gauge field produced by the "external" current j^. The calculations are lengthy 
and their details are adduced in Appendix 3. Here, we present only the final answers. 

The components of the longitudinal propagator A[,^''(t2, r, ry) are already obtained in the form of the three- 
dimensional integrals ( 2.43| ). A*^-^) depends on the differences of the curvilinear spatial coordinates, r = ri — r-i 
and r] = r]i — r]2, but not on the difference of the temporal arguments ti of the field and T2 of the source. Intro- 
ducing the shorthand notation for the distance in the (a;y)-plane, r± — and for the full distance i?2 = R(j'2) = 
[(ri — r2)^ -|- t| sinh^(77i — ?72)]^^^ between the two points of the surface T2 ~ const, we obtain: 



A(i) = - 



0{t2 



92 



47r 



d-qdx' 
. 0{t2 



■(i2 C0th|?7|) 



4:71 dx'^dx^ 
9{t2 - r^) 



1 



L2 coth — In 
92 / 



47r cosh 77 977 \ sinh 1 77 1 



^ 7] coth 77 — 1 
sinh^ 77 



vi 

2 



rj_ cosh 7] 
sinh^ 1 7; I 



T2-R2 cosh 7/ 
sinh^ 77 



(4.1) 



where L2 — L{t2) — ln[(T2 sinh |7;| + R2)/r 
A(^) = 



^2 



A 



An 



After the derivatives are evaluated, most of the logarithms here vanish: 



T2 cosh 77 



R2 



2x'^x^ T2 cosh 77 



Rl 



(t2 - rj^) T2 sinh 77 t| 



All 

B{t2- 



An 



R2 
7] coth 77 



Rj 



6{t2 — r_i_) x^ rl sinh 77 



An 



Rl 



1 T2 cosh 77 



sinh 77 



i?2 sinh 77 



3- ^ 

i?9 



2 cosh 77 
sinh"^ 1 77 1 



(4.2) 



By examination of Eq. (2.42), one may see that after replacement of T2 by n the same kernel, A[^' (ti, r*, 77), determines 

the components Em^ (ji ) of the longitudinal part of the electric field via the components (n ) of the current at the 

same time. 

Similar calculations have led to the expression (2.20) for the electric field of the static Coulomb source. The kernel 
Km of the instantaneous potential has the following components, 



'^r = cosh 7;— 



(ti — r±) x^ T2 cosh 77 sinh 1 7/ 1 



An 



IC,, 



0{ti 



An 



Rl 

tanh 1 77 1 



Rl sinh 1 77 1 



sinh^ 77 



Li 



(4.3) 



where Ri — R{ti) and Li = L{ti). These propagators do not respect the light cone, but have a remarkable property 
that the longitudinal fields at the surface of the constant proper time r do not exist at the distances r± from their 
sources that exceed t. This establish the upper limit for the possible dynamical correlations between the longitudinal 
fields in the {xy)-p\a,ne. 



V. GLUON CORRELATORS IN THE CENTRAL RAPIDITY REGION AND NEAR THE LIGHT 

WEDGE 

Our next step is to compare the correlators of the gauge = with the similar correlators in three other gauges, 
A" = 0, A'^ = and — 0. We shall start with the simplest on-mass-shell Wightman function A^q . This type 
correlators, Aq^, Aq'^ and A^'' share the same polarization sum of the free gauge field. They correspond to the 
densities of the final states of the radiation field and are important for various calculations. The same polarization 
sum appears in expressions for the transverse part of the propagators, A^^^, A(^^^, Agg and Ai^. For our immediate 
purpose we shall include the projector d'"' of the gauge A'^ = Q into the formal Fourier representation. 
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/(Pk 



-ik{xi—X2) 



(5.1) 



with the "extraneous" dependence of the Fourier transform on the time and spacial coordinates. This dependence 
disappear in some important hmits. Therefore, we discover the domains where the wedge dynamic simphfies and 
describe the processes which are approximately homogeneous in space and time. These domains are: (i) the central 
rapidity region, r^i 2 «C 1 (or 2 ~ 0)i where the projector in the integrand of Eq. (5.1) is 



ku 



kt^k" 
{kuf 



(5.2) 



with the gauge-fixing vector, = (1, 0, 0, 0), which approximately coincide with the local normal to the hypersurface 
T — const] and (ii), the vicinities of two null-planes, 77 — > ±00 (or 0) where 



d^''{n±,k)^~g'^'' + 



k^'n^J^ + k'^n'^ 
{kn±) 



with the null-plane vectors — (1, 0, 0, =Fl) . 

Eqs. (3.(:)-( [3j| ) almost fit our needs. In all three 0, fc°a;° > 1, as well as x 



0, fc^x 3> 1, the functions /i and /2 can be approximated by the following expressions 



/i « i tanh(6' - ,^)e-*fc_LTcosh(e-r,)+*fcr 



-ikx 



k 



-ikx 



k+3 



(5.3) 

0, k~x+ > 1 and 
(5.4) 



^ — ik±T cosh(9—ri)-\-ikr ^—ikx ^ — ikx 
f2 ~ ik±T — — = ik'^i t'^ , n n ,11 — T^— ; . (5.5) 

cosh(6'-?7) ^ k^x"-k^x^ ^ k+x-+k-x+ ^ ' 

(We have omitted the time independent terms in /i and ji which set the potentials of the mode u^^^ to zero at r = 0. 
This kind of terms would correspond to the residual gauge symmetry and is not kept in the axial and the null-plane 
gauges as well.) 

Transformation of the correlator A'™(1, 2) to the Minkowsky coordinates is carried out according to the formula, 

D^^'^ix^.x^) = af (xi)g''(a:i)A,„,(ui,M2)g""(a:2)a5;(x2) (5.6) 
where the matrix of the transformation is defined in a standard way, 

at{x) = —, a'l{x)^x\ al{x)^x\ < = . (5.7) 

These are the only components of the tensor a'' (a;) which participate the transformation. In this way we obtain 

i:»°"(l, 2) = 2;?a;^A'"'(l, 2); D'^'^il, 2) = a;j2;°A'"'(l, 2) 
^30(1, 2) = xlxl^^^{l, 2); D^^l, 2) = x\xl^^^{l, 2), 
Z?"'-(l, 2) = x?A'"-(l, 2); D^\l, 2) - a;^A'^''(l, 2); 
D^'-il, 2) x?A'"-(l, 2); £''■3(1, 2) ^ a;^A'-''(l, 2), 

£''■^(1,2) = A'-^(l,2) . (5.8) 

Every additional factor g'''' = r^^ finds a counterpart which prevents a singular behavior at t = 0. In the above 
approximation, the expression for the A''''(xi, 2:2) component of the correlator has the form: 

r d^k p^^-rk(x^~X2) ^ r ^J^2^^-^k(x^-X2) 

[Xi,X2) - J (2^)32fc0 (p^O _ fc32.3)(p^0 _ ^^.3) - J (27r)32fc0 (fc + a;- + k-x+){k+x^ + k-x+) ' 

Therefore, in the limit of 2 — > we obtain that 0, while d^^{k,u) k'j^/kQ, thus reproducing the 

corresponding components of the gauge A'^ ~ 0. The other components are reproduced one by one as well, and one 
can expect smooth transition between the gauge of the the wedge dynamic and the local temporal axial gauge of the 
reference frame co-moving with the dense quark-gluon matter created in the collision. 
In the limits of x^2 ^ we obtain that 
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£,03^ £,30^ £,33 Ji_ a , (5.10) 

and 

DOO, i)30^ ^33 _ = ^ , if ^ . (5.11) 

These limits, after they are found for all components, lead to the well known expressions of the projectors d'^'^{n, k) 
in the null-plane gauge: the gauge = in the vicinity of a;+ = 0, and A~ = in the vicinity of x~ = 0. Therefore 
we obtained an expected result: in the limit of the light-cone propagation the gauge = recovers the null-plane 
gauges = and = 0. 

Some remarks are in order. First, the conception of the structure functions relies heavily on the null-plane dynamics 

which essentially uses these gauges. For two hadrons (or two nuclei) we have two different null-plane dynamics which 
do not share the same Hilbert space of states. Now wc have an important opportunity to describe both nuclei and 
the fields produced in their interaction within the same dynamic and the same Hilbert space. Second, one may trace 
back the origin of the poles {ku)~^ in the polarization sums of axial gauges (uA) = and see that they appeared in 
course of approximation of the less-singular factor [A;_l cosh(0 — 77)] ~^ in various limits of the propagator of the gauge 
A'' = 0. 

Further, contrary to a naive expectation to obtain the gauge A'^ = at a;~ =0 and the gauge A~ = at a;+ = 0, we 
obtained them in the opposite correspondence. First of all, let us notice that the result is mathematically consistent. 
Indeed, the gauge condition A'^ = may be rewritten in the form, 

= ^(A+e-^ + A-e") =0 . (5.12) 

Thus the limit of 77 —+ 00 {x~ — * 0) indeed leads to A~ = and the limit oi rj ^ —00 (a:+ 0) leads to A'^ = as 
the limit gauge conditions. Recalling that 

A" = ^{A+e-'' -A-e'') = A+e-'' = -A-e'> , (5.13) 

we immediately realize that in the vicinities of both null-planes the tangent component A^ = 0. This fact has a very 
simple geometrical explanation: the normal and tangent vectors of the null plane arc degenerate. Once A'^ = 0, we 
have A^ = and consequently, A^ = and A~ = at r/ ^ ±00. This result naturally follows from the geometry of 
the system of the surfaces where we define the field states. These are subject for dynamical evolution in the direction 
which is normal to the hypersurface. 

The perturbative QCD-evolution for the object which moves in the a::+-direction and experience the deep inelastic 
scattering, begins only after it has passed the collision point a;+ = = 0. The so-called pre-collision dynamics cannot 
start earlier: the colliding hadron and electron are Lorentz contracted to a small longitudinal size. It is important 
that the gauge field correlators of the null-plane gauges were obtained as the sums over the states which emerge from 
the collision point. At r ^ the stationary phases of all in- wedge plane- wave modes are uniformly spread along both 
null-planes. If no interaction happens, then the decomposition remains virtual. The partial waves will keep their 
coherence and assemble into the incoming hadron and the electron travelling in the initial directions. It takes finite 
time to distort the phase balance between the partial waves of the initial states. Before the finite distortion develops, 
the wave packet of the hadron partially goes through the assembling towards the state of normal localized hadron. 
Only a fraction of the hadron interacts with the structureless electron. Therefore, the perturbative QCD-evolution of 
the hadronic wave function indeed takes place in the vicinity of the null-plane a;"*" = 0, where the limit of the gauge 
A'^ = is A~^ = 0. It is exactly the domain of large negative r] where the interaction which excites the state of the 
most hard back scattered quark is localized. This qualitative picture is in compliance with the well known fact that 
the QCD evolution equations are adequate only when x is not too close to 1. 

It is expedient to emphasize similarity of our approach to the classical spectral analysis of the transient processes 
and to the standard scheme of the quantum-mechanical measurement. Thoiigh wc decompose the wave packet of 
the incoming hadron in terms of the states which are immediately before and after the collision are spread over the 
entire surface of the light wedge, the phases and amplitudes of this decomposition are balanced in such a way that 
before the collision the whole packet is strongly localized in the vicinity of the null-plane x~ = 0. One needs a strong 
interaction in the vertex in order to break down this tiny balance and filter out an unusual field configuration. 
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APPENDIX 1. Modes of the free gauge field 

Here, we shall obtain the complete set of the on e-pa rticle solutions to the homogeneous system of the Maxwell 
equations with the gauge A'^ = 0, that is Eqs. ( ^.lC| )-( ^^ ). This gauge condition expHcitly depends on the coordinates, 
thus introducing effective non-locality in the path integral that represents the action. Therefore it becomes impossible 
to invert the differential operators using the standard symbolic methods. The knowledge of the one-particle solutions 
becomes necessary in order to find the Wightman functions of the free vector field, to establish the the explicit form 
of the field commutators, and to separate the propagators of the transverse and the longitudinal fields. It is natural 
to look for the solution in the form of the Fourier transform with respect to the spacial coordinates: 

A,{x)^ f dv [ d^k e'^'^e'^^ A,{k,v,T) (Al.l) 



Then the system of the second order ordinary differential equations for the Fourier transforms becomes as follows: 

ly^ —* —' lyk —* 

[Tdl+dr + — +Te,Mk,U,T)-Tk.^kyAy{k,V,T) Ar^{k , U,t) = G , (Al.2) 

— * i/^ — * jyk ~* 

- Tk^kyA^{k,i',T) + [Td'^ + dr ^ \- Tkl]Ay{k,iy,T) -AJk,i',T) = , (Al.3) 

r T 



-'^A^{k,,y,T)-^Ay{k,,,,T) + [-d^^-^dr + -kl]A^{k,,y,T)=0 , (A1.4) 
In this form the system is manifestly symmetric and self-adjoint. An additional equation of the constraint reads as 



C{k,iy,T) — —vdrA,, + Tdr[kxAx{k,i',T) + kyAy{k,i^,T)] — . (Al.5) 

T 

Let us rewrite the homogeneous system of the Maxwell equations in terms of the variables 

^ ^ dxAa; + dyAy, ^ dyA^, - dxAy, andA==A,,. (Al.6) 

One immediately sees that the equation for the Fourier component ^'(fc, ly, t) of the longitudinal magnetic field ^>{x) 
decouple: 

[dl + -dr + ^ + kl]-^{k,i^,T) = (A1.7) 

T T 

Then the other two equations of motion take shape 

[T'd^, + Tdr + i^']'i>f:^jT)-zi.klAj:jT) = , (A1.8) 

[d'^-^d, + kl]AgjT)+ti.<i>j:jT) = , (A1.9) 
The additional constraint equation can be conveniently rewritten as 

C{k, V, t) = ^a,A,~^(r) + Tdr^j:jT) = . (ALIO) 
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This is an independent equation. However, the conservation of the constraint along the Hamiltonian time t is 
a consequence of the equations of motion, and it can be employed to obtain the independent equations for the 
components of the vector field. This is easily done in terms of the auxiliary functions. 



nJ^)=^^kJ^) and aj:jT)=T ^Aj:At), 



(Al.ll) 



which are directly connected to the "physical" components of the electric field, 5™ = -y/— gg™'yli 



(A1.12) 



As a result, we obtain that the functions (pj: ^(r) and ^(r) obey inhomogeneous Bessel equations. 



(A1.13) 



(A1.14) 



(A1.15) 



where Cip and Ca are arbitrary constants. Now we may cast the solution of these equations in the form of the sum of 
the partial solution of the inhomogeneous equation and a general solution of the homogeneous equation. 



(A1.16) 



*rr(l) 



(A1.17) 



(1) 



(A1.18) 



where 3^.1/(2:) is the so-called Lommel function 

Furthermore, it is useful to notice that the system of the Maxwell equations ( ^.10| )-( p.ll| ) also has an infinite set of 
the T-independent solutions of the form 



(A1.19) 



where x is an arbitrary function of the spacial coordinates 77 and r. Thus, they are the pure gauge solutions of the 
Abelian theory, compatible with the gauge condition. 



In order to find the coefficients one should integrate Eqs. (Al.17) and (Al.l^) with respect to the Hamiltonian 



time T, thus finding the functions $ and A. Next, it is necessary to solve Eqs.(Al.6) for the Fourier components of 
the vector potential and to substitute them into the original system of Eqs. ( Al.2)-(A1.5). Using functional relations 



from Appendix 2, one obtains that c -I- 1^7 = and Ca — vci, 







One of the solutions, (already normahzed according to Eq.(p.9|)) is found immediately: 



-k^ I i/i'U*i^)e''"'^'^' 




(A1.20) 



Initially, the components of the vector mode V^"^^ appear in the following form required by the convergence of the 
integral, 



Rx\iv^kx_T\s) - iz^[e'^''/Vsinh(7riy/2)] 



,ykrR%^.,,{k^T\S) 



However, it can be gauge transformed to the more compact form. 
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The third solution, the last one by the count of the nonvanishing components of the vector potential in the gauge 
= has the following form, 

yi^)/^.) ^ ( krQ-l,^Ak±T) \ ^,,^+,kp (A1.22) 

The modes V'^^^ and F^^) are the normalized solutions of the Maxwell equations. They are orthogonal and obey 
the normalization condition, 

(V^f\v^:^)-S(.-.')6(k^k'), iy^,V^ = . (A1.23) 



which can be easily verified by means of Eq. (A2.5). The norm of these solutions is given by the Eq. ( |2.9D . A 
normalization coefficient of the mode F*-^-* is not defined as this mode has a zero norm. It is also orthogonal to y*^^^ 
and 

(Vi^), ) yi^') = = . (A1.24) 

k,v k' ,u' k.u k,v k.u k,v 

Thus this mode drops out from the decomposition of the free gauge field. 

The conservation of the constraint can be obtained as a consequence of the Eqs. ( A1.12| ) and (Al.13) in the form, 

^dr[VuA^) + ^ ^drCj:u{T) = , (A1.25) 

which reassures us in consistency between the dynamic equations and conservation of th e Gaus s law constraint. 

One can explicitly check that the modes V^^'' and F^^-' obey the constraint equation ( ALlOj ), which expresses the 
Gauss law. The mode V^^'^ does not. This mode corresponds to the longitudinal field which cannot exist without the 
source. 

APPENDIX 2. Mathematical miscellany 

This appendix contains a list of mathematical formulae for the functions which appear in various calculations in 
the body of the paper and Appendix 1. The components of the vector field are expressed via two types of integrals. 
The first of them was studied in Ref. |^ : 

= Jx''Hl^\x)dx = .t[(^ + - l)i/(^')(x)S^_i,,_i(:r) - hI'1,{x)S^,,{x)] , (A2.1) 

where S^,,^ stands for any of the two Lommel functions, s^^i, or S^^i, [ Whenever we omit the indicator |S), the 

function i?jf^]y(a;|s) is assumed.] The second type of integrals, 

Q^,,.ix) = / (A2.2) 
Jo 

is a new one. The functions i?jf^t(a;|S) are introduced as the indefinite integrals. The preliminary choice of the 
lower limit and, consequently, between s^^^ and 6*^1,1/, is motivated by the requirement of the convergence and regular 
behavior. One can easily prove that 

(2\ (2\ -p- 7ri//2 

i;y.,„WS)-«(i.„(.|.).-i^. 

R\%,MS) R\%.AAs) = . (A2.3) 

smh(7rz^/2j 

We often use the following relation between the Lommel functions 

Si,..(fciT) = 1 - v^S-i^Uk±T) . (A2.4) 
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From the integral representations (A2.1) and (A2.2). it is straightforward to derive the functional relations 



R 



U) 



-l.iv 



(A2.5) 



The Wronskian of the Hankel and Lommcl functions, 



(A2.6) 



(A2.7) 



is necessary to prove orthogonality of V^*^^' and y^^^. To prove ( A2.7 ), one should use the following representation for 
the Lomniel function, 



?(2 



(2) /^N 7,(1) 



Ai 



which follows from Eq. (A2.1) and its consequence. 



(A2., 



(A2.9) 



In order to prove relation ( 2.33 ) one should use representation ( A2.1 ) for the functions Ri^J, and the Wronskian of 
two independent Hankel functions. The proof of relation (2.44) begins with replacing the functions _R^| by R^il^ 
by means of Eq. ( |A2.5| ). The final result follows from Eq. ( A2.9D and ( |A2.6| ). 

APPENDIX 3. Calculation of the longitudinal part of the propagator 



The kernels (2.45) and (2.47) of the longitudinal and instantaneous parts of the propagator are given in the form of 
the three-dimensional Fourier integrals dvcPk. Here, we describe major steps of calculations which lead to Eqs. (4.1 
and (O). 



We permanently use the following integral representation for the Hankel functions, 

^ J —oo 



(A3.1) 



which allows one to calculate many integrals by changing the order of integration. The Lommcl function Si^ii, has a 
similar representation, 



poo 

Si,iv{x) = X / coshucosi/M e"^™''"du , 

^0 



Integrating it by parts, and using Eq. ( A2.4 ), we find the integral representation for S-i^^, 

/>oo 



(A3.2) 



(A3.3) 



We start with integral representation (A2.2) of the functions Q±i^ii, and perform integration over v. To compute the 
integrals from the function si it can be conveniently decomposed in the following way. 



2 sinh(7r!^/2) ' 



(A3.4) 



which allows one to find 



Si^^^{k^T)e^'"^dv = 7rfc_Lrcosh?7e^'^--L^^'"''l"l , 

/OC 
vS^x,^,{k^T)e"'''dv = msignTye-^'^"*^'"^^!"! , 
-OO 



(A3.5) 
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The similar Fourier integrals from the function hi^ are calculated using the representation (A3.1) for the Hankel 
functions and and the integral, 



7rj//2 



1 



sinh(7rz^/2) 



de- 



cosh 



— ^ ^ sin [fcj_r cosh (6* — ri)] 
cosh^e ^ ^ 



This yields, for example, 

/oo 
di^e'-''^h,,{k^T) 
-oo J-oo cosn a 

After integration over v we obtain the following integral for the Alt'* 



(2^~fc2-^ 



— TrfciTCOshne-'^^^'^'"'^!''! 



2 — sin[fcxT cosh(0 — t])] 

cosh 6* 



(A3.6) 



(A3.7) 



(A3.. 



and similar integrals for the other components. The first term in this formula is calculated in the following way. After 
integration over r we continue: 



A' = 



' 87r2 



coth \7j\ I 4^ e'^^ll - e 



'k±T2 sinh \ ri\' 



47r 



coth 1 77 1 



^Jo(fc±rx)[l-e-^^"^^'"''l"l] = -Mlcoth|7/|ln 
k± Ait 



T2 sinh 1 77! 



'1 + '''2 sinh^ 77 



(A3.9) 



To work out he second term, one should introduce kz — fcj_sinh0 and fco = fc^coshf? = |k| and join d^kd9 in one 
three-dimensional integration d^k. With t = r cosh 77, r = (x, y, t sinh?]), this leads to 



A" = 



drd, n dr f d^k 



— / — ^e''^'"sinfcnt = 



drds 



7,3 



47r 



T cosh 77 



+ r2 sinh^ 



drds 
An 



0(r± ~ T2) cothj?;! In 



T2 sinh 1 7; I 



t| sinh 77 



+ e{T2-ri_) In 



T2 



(A3. 10) 



Adding (A3. 9) and (A3. 10), we obtain the first of the equations (fl.l[). 
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